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INTRODUCTION 


This  report  represents  the  final  report  for  the  research 


carried  out 


under 


ARPA-ONR  Contract 


Number 


N00014-75-C-0951-NR  049-328.  ^During  the  course  of  this  contract, 
there  have  been  several  major  task  areas: 

1.  The  development  of  homomorphic  signal  processing  tech¬ 
niques  and  their  application  to  the  development  of  a  homomorphic 
vocoder  and  other  signal  processing  applications. 

2.  The  development  and  implementation  of  techniques  for 
enhancement  and  bandwidth  compression  of  degraded  speech. 

3.  The  development  and  evaluation  of  techniques  for  pro¬ 
cessing  of  multidimensional  signals  and  the  application  of  these 
techniques  to  Image  processing,  and  the  processing  of  other 
multidimensional  data  sets.  — 

4.  The  development  of  techniques  and  the  Implementation  and 
evaluation  of  systems  for  speed  transformations  of  speech. 

\ 

In  addition,  in  the  course  of  this  contract,  there  have'  been 
a  number  of  other  studies  carried  out  relating  to  the  development 
of  signal  processing  techniques,  as  well  as  their  implementation 
and  application  to  particular  problems  of  interest. 

Throughout  the  course  of  this  contract,  we  have  reported  the 
details  of  our  various  research  efforts  in  the  form  of  Journal 


V 
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articles  and  In  some  cases  technical  reports.  At  the  end  of  this 
final  report  Is  a  complete  listing  of  the  technical  reports  and 
Journal  articles  reporting  work  carried  out  under  the  contract. 
A  number  of  the  papers  supported  under  this  contract  and  listed 
In  the  bibliography  received  major  IEEE  awards.  In  the  following 
sections,  we  give  a  brief  overview  of  the  work  carried  out  under 


the  various  categories  indicated  above. 


i 
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HOMOMORPHIC  SIGNAL  PROCESSING 

A  major  focus  of  the  work  under  this  contract  has  been  the 
development  of  homomorphic  signal  processing  techniques  directed 
principally  at  the  development  of  a  low-cost  homomorphic  vocoder. 
As  part  of  this  work,  a  number  of  basic  new  algorithms  were  deve¬ 
loped  for  homomorphic  signal  processing  including  the  development 
of  a  new  phase  unwrapping  algorithm  and  the  development  of  tech¬ 
niques  for  homomorphic  processing  using  both  the  conventional  and 
sliding  chirp  z-transform.  Also  out  of  this  study  emerged  a 
synthesis  of  the  techniques  of  homomorphic  filtering  and  linear 
prediction  analysis.  This  new  modified  technique  has  been 
labeled  as  homomorphic  prediction. 

Much  of  the  emphases  on  the  application  of  homomorphic  tech¬ 
niques  to  speech  compression  has  centered  around  the  incor¬ 
poration  of  phase.  To  accomplish  this,  a  set  of  detailed  tech¬ 
niques  for  phase  estimation  with  particular  application  to  speech 
analysis-synthesis  was  explored  and  a  number  of  iterative  tech¬ 
niques  were  developed  to  accomplish  the  phase  estimation.  The 
importance  of  phase  in  relation  to  speech  processing,  as  well  as 
the  processing  of  other  signals,  was  a  central  theme  throughout 
this  work.  In  addition  to  the  development  of  a  system  for  the 
incorporation  of  phase,  a  more  general  theory  relating  to  the 
importance  of  phase  in  both  speech  and  pictures  was  developed. 
This  basic  and  fundamental  work  was  reported  in  an  invited  IEEE 
paper  published  in  May  1981. 


/ 
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Another  modified  homomorphic  technique  which  was  developed 
and  explored  in  considerable  detail  is  referred  to  as  spectral 
root  homomorphic  deconvolution.  In  addition  to  the  development 
of  the  theory,  the  spectral  root  homomorphic  vocoder  was  imple¬ 
mented  and  evaluated. 

In  addition  to  its  application  to  speech  processing,  work  was 
carried  out  to  explore  and  develop  the  application  of  homomorphic 
filtering  in  the  context  of  seismic  data  processing. 

SPEECH  ENHANCEMENT  AND  BANDWIDTH  COMPRESSION 

A  second  major  project  during  the  course  of  this  contract  was 
the  development,  implementation  and  evaluation  of  techniques  for 
enhancement  and  bandwidth  compression  of  degraded  speech.  The 
principal  accomplishments  on  this  project  were  reported  in  detail 
in  an  Invited  review  paper  published  in  the  IEEE  Proceedings  in 
December  1979*  The  system  which  was  developed  and  implemented 
centered  around  the  all-pole  modeling  of  speech  taking  into 
account  the  noisy  observations.  In  addition  to  the  implemented 
system  being  better  in  performance  than  other  previous  systems, 
a  considerable  body  of  important  and  useful  theory  relating  to 
the  modeling  of  signals  which  have  been  degraded  by  noise  was 
developed. 
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MULTIDIMENSIONAL  SIGNAL  PROCESSING 

Our  work  on  multidimensional  signal  processing  had  a  very 
strong  component  relating  to  the  development  of  the  theory  for 
multidimensional  filtering  and  spectral  analysis.  A  number  of 
new  multidimensional  filter  design  techniques  were  developed  and 
are  now  widely  used  throughout  the  community.  Much  of  our  early 
work  in  this  area  was  published  in  an  invited  paper  in  the  IEEE 
Proceedings  on  two-dimensional  digital  filtering.  Important 
results  were  also  obtained  under  this  project  on  multidimensional 
spectral  analysis.  A  basic  new  algorithm  for  two-dimensional 
maximum  entropy  power  spectrum  estimation  was  developed  which 
resolves  many  of  the  algorithm  problems  previously 
encountered  in  relation  to  multidimensional  power  spectrum  analy- 
sis.  A  new  class  of  algorithms  was  also  developed  for  the 
numerical  evaluation  of  the  hankel  transform  which  corresponds  to 
the  two-dimensional  Fourier  transform  of  a  circularly  symmetric 
function.  This  latte  algorithm  was  motivated  in  part  by  a 
problem  relating  to  the  measurment  of  the  plane  wave  reflection 
coefficient  of  the  ocean  bottom  and  is  now  generally  considered 
to  be  one  of  the  most  efficient  algorithms  for  the  computation  of 
the  hankel  transform. 

The  theory  that  has  been  developed  under  this  project  for 
multidimensional  signal  processing  has  been  applied  in  a  number 
of  areas.  One,  as  indicated  above,  is  in  problems  relating  to 
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ocean  acoustics.  Also,  much  of  the  work  has  been  implemented  in 
relation  to  problems  of  image  enhancement  and  restoration 
including  the  restoration  of  speckle  images. 

SPEED  TRANSFORMATIONS  OF  SPEECH 

During  this  contract,  we  developed,  implemented,  and  investi¬ 
gated  in  detail  techniques  for  speed  transformations  of 
speech,  in  particular,  techniques  for  speeding  up  and  slowing 
down  speech  while  retaining  very  high  quality.  In  addition  to  a 
successful  implementation  which  resulted  in  a  system  of  con¬ 
siderably  higher  quality  than  had  ever  previously  been  achieved., 
a  basic  new  theory  resulted  on  the  short  time  Fourier  analysis  of 
speech  and  the  use  of  this  theory  in  the  implementation  of  the 
digital  phase  vocoder.  r 

GENERAL  SIGNAL  PROCESSING 

In  part  as  an  outgrowth  of  the  projects  described  above  a 
number  of  other  important  studies  were  carried  out  under  this 
contract.  A  theory  was  developed  for  the  analysis  of  linear 
digital  networks  which  was  published  in  an  IEEE  invited  paper. 
Work  on  computer  architectures  for  signal  processing,  also 
published  aq  an  IEEE  Proceedings  invited  paper,  was  carried  out. 
More  recently  work  under  this  contract,  motivated  in  large  part 
by  our  previous  work  on  the  importance  of  incorporating  phase 


v 
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information  in  vocoders,  was  carried  out  on  the  importance  of 
phase  in  signals.  Out  of  this  study  has  emerged  a  rich  and 
potentially  very  significant  theory  on  signal  reconstruction  from 
either  phase  or  magnitude  information  in  the  time  domain  or  fre¬ 
quency  domain. 


v 
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PUBLICATIONS  1975  -  1981 
under  Contract  N00014-75-C-0951-NR  049-328 

JOURNAL  ARTICLES,  REPORTS  AND  BOOKS 


1975 


1.  D.  Dudgeon,  "Existence  of  Cepstra  for  Two-Dimensional 
Rational  Polynomials,”  IEEE  Trans,  on  Acoustics,  Speech  and 
Signal  Processing,  vol.  ASSP-23,  no.  2,  pp.  242-243,  April  1975 • 

2.  R.  Crochiere  and  A.  Oppenheim,  "Analysis  of  Linear  Digital 
Networks,"  IEEE  Proceedings  Special  Issue  on  Digital  Signal 
Processing,  invited,  vol.  63,  no.  4,  pp.  581-595,  April  1975. 

3.  R.  Mersereau  and  D.  Dudgeon,  "Two-Dimensional  Digital 
Filtering,"  IEEE  Proceedings  Special  Issue  on  Digital  Signal 
Processing ,  invited,  vol.  63,  no".  4,  pp.  610-623,  April  1975* 

4.  J.  Allen,  "Computer  Architecture  for  Signal  Processing," 
IEEE  Proceedings  Special  Issue  on  Digital  Signal  Processing, 
invited,  vol.  63,  no.  4,  pp.  624-633,  April  1975. 

5.  D.  Dudgeon,  "Two-Dimensional  Recursive  Filter  Design  Using 
Differential  Correction,"  IEEE  Trans,  on  Acoustics,  Speech  and 
Signal  Processing,  vol.  ASSP-23,  no.  3,  pp.  264-267,  June  1975. 

6.  R.  H.  Frazier,  S.'  Samsam,  L.  D.  Braida  and  A.  V.  Oppenheim, 
"Enhancement  of  Speech  by  Adaptive  Filtering,"  RLE  Progress 
Report  No.  116,  pp.  177-181,  July  15,  1975. 

7.  E.  Singer  and  A.  Oppenheim,  "Amplitude  Filtering  Using  a 
Time  Variable  Transformation,"  RLE  Progress  Report  No.  116,  pp. 
181-188,  July  15,  1975- 


1976 

8.  A.  V.  Oppenheim,  M.  R.  Portnoff,  "Speed  Transf ormations  of 
Speech,"  RLE  Progress  Report  No.  117.  p.  75,  January  1976. 

9.  A.  Oppenheim  and  E.  Singer,  "Enhancement  of  Lowpass  Filtered 
Speech,"  RLE  Progress  Report  No.  117,  pp.  75-76,  January  1976. 

10.  A.  P.  Holt,  W.F.G.  Mecklenbrauker  and  A.  V.  Oppenheim, 
"Design  and  Implementation  of  Variable  Cutoff  Digital  Filters," 
RLE  Progress  Report  No.  117,  p.  76,  January  1976. 

11.  Q.  E.  Kopec,  A.  V.  Oppenheim,  "Speech  Analysis  by 
Homomorphic  Prediction,"  RLE  Progress  Report  No.  117.  pp.  76-77, 
January  1976. 
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12.  A.  V.  Oppenheim,  J.  M.  Tribolet,  '’Applications  of 
Homomorphic  Filtering  to  Seismic  Data  Processing,”  RLE  Progress 
Report  No.  117.  p.  77,  January  1976. 


13.  H.  W.  Schuessler,  "A  Stability  Theorem  for  Discrete 
Systems,"  IEEE  Trans,  on  Acoustics,  Speech  and  Signal  Processing, 
vol.  ASSP-24 ,  no.  1,  pp.  87-89,  February  1976. 

14.  A.  V.  Oppenheim,  W.  F.  G.  Mecklenbrauker  and  R.  M. 
Mersereau,  "Variable  Cutoff  Linear  Phase  Digital  Filters,"  IEEE 
Trans,  on  Circuits  and  Systems ,  vol.  CAS-23,  no.  4,  pp.  199-203, 
April  197(>* 

15.  M.  R.  Portnoff,  "Implementation  of  the  Digital  Phase  Vocoder 
Using  the  Fast  Fourier  Transform,"  IEEE  Trans,  on  Acoustics, 
Speech  and  Signal  Processing,  vol.  ASSP-24 ,  no.  3,  pp^  243-248 , 
JuneTgfST 

16.  R.  M.  Mersereau,  W.  F.  G.  Mecklenbrauker  and  T.  F.  Quatleri, 
Jr.,  "McClellan  Transformations  for  Two-Dimensional  Digital 
Filtering:  I.  Design,"  IEEE  Trans,  on  Circuits  and  Systems, 
vol.  CAS-23,  no.  7,  pp.  405-414,  July  19T57 

17.  W.  F.  G.  Mecklenbrauker  and  R.  M.  Mersereau,  "McClellan 
Transformations  for  Two-Dimensional  Digital  Filtering:  II. 
Implementation,"  IEEE  ,  Trans .  on  Circuits  and  Systems,  vol. 
CAS-23,  no.  7,  pp.  414-422,  July  ISlV. 

18.  A.  V.  Oppenheim,  G.  Kopec  and  J.  Tribolet,  "Signal  Analysis 
by  Homomorphic  Prediction,"  IEEE  Trans,  on  Acoustics,  Speech  and 
Signal  Processing,  vol.  ASSP-24,  no~.  4 ,  pp.  327-332,  August  1976. 

19.  A.  P.  Holt,  "A  Variable  Cutoff  Filter  Using  Digital 
Frequency  Warping,"  S.M.  Thesis,  M.I.T.,  Dept,  of  Elec.  Eng.  and 
Comp.  Science,  September  1976 . 


1977 

20.  D.  B.  Harris,  J.  H.  McClellan,  "Reconstruction  of  Velocity 
Structures  from  Teleseismic  First  Arrival  Times,"  RLE  Progress 
Report  No.  119,  p.  98,  January  1977. 

21.  J.  M.  Tribolet,  A.  V.  Oppenheim,  "Application  of  Homomorphic 
Filtering  to  Seismic  Data  Processing,"  RLE  Progress  Report  No. 
119 .  P*  99,  January  1977. 

22.  J.  S.  Lim  and  A.  V.  Oppenheim,  "Enhancement  of  Degraded 
Speech,"  RLE  Progress  Report  No.  119.  pp.  99-100,  January  1977. 
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23  •  E.  Singer  and  A.  V.  Oppenheim,  "Enhancement  of  Lowpass 
Filtered  Speech,”  RLE  Progress  Report  No.  119,  pp.  100-101, 
January  1977 •  . . 

24.  T.  F.  Quatieri  and  A.  V.  Oppenheim,  "Speech 
Analysis-Synthesis  Based  on  Homomorphic  Filtering  and  CCD 
Technology,"  RLE  Progress  Report  No.  119.  p.  101,  January  1977. 

25.  M.  R.  Portnoff  and  A.  V.  Oppenheim,  "Speed  Transformations 
of  Speech  Signals,"  RLE  Progress  Report  No.  119,  pp.  101-102, 
January  1977.  '  . . . 

26.  G.  E.  Kopec,  A.  V.  Oppenheim  and  J.  M.  Tribolet,  "Speech 
Analysis  by  Homomorphic  Prediction,"  IEEE  Trans,  on  Acoustics, 
Speech  and  Signal  Processing,  vol .  ASSP-25,  no.  1,  pp.  40-49, 
February  1977. 

27-  J.  M.  Tribolet,  "A  New  Phase  Unwrapping  Algorithm,"  IEEE 
Trans .  on  Acoustics,  Speech  and  Signal  Processing,  vol.  ASSP-25, 
no.  2,  pp.  170-177,  April  1977. 

28.  J.  M.  Tribolet,  "Seismic  Applications  of  Homomorphic  Signal 
Processing,"  Sc.D.  Thesis,  M.I.T. ,  Dept,  of  Elec.  Eng.  and  Comp. 
Science,  May  1977. 

29.  D.  B.  Harris  and  R.  M.  Mersereau,  "A  Comparison  of 
Algorithms  for  Minimax  Design  of  Two-Dimensional  Linear  Phase  FIR 
Digital  Filters,"  IEEE  Trans,  on  Acoustics,  Speech  and  Signal 
Processing,  vol.  ASSP-25,  no.  6.  pp.  492-500,  December  1977. 


1978 

30.  J.  M.  Tribolet,  "Applications  of  Short-Time  Homomorphic 
Signal  Analysis  to  Seismic  Wavelet  Estimation,"  Geoexploration , 
vol.  16,  pp.  75-96,  January  1978. 

31.  A.  V.  Oppenheim,  G.  V.  Frisk  and  D.  R.  Martinez,  "An 
Algorithm  for  the  Numerical  Evaluation  of  the  Hankel  Transform," 
IEEE  Proceedings  Letters,  vol.  66,  no.  2,  pp.  264-265,  February 
197F. 

32.  J.  S.  Lim  and  A.  V.  Oppenheim,  "All  Pole  Modeling  of 
Degraded  Speech,"  IEEE  Trans,  on  Acoustics,  Speech  and  Signal 
Processing,  vol.  ASSP-26,  no,  3,  pp.  197-210,  June  197 BT 

33*  J.  M.  Tribolet,  "Applications  of  Short-Time  Homomorphic 
Signal  Analysis  to  Seismic  Wavelet  Estimation,"  IEEE  Trans,  on 
Acoustics ,  Speech  and  Signal  Processing,  vol.  ASSP-26,  no.  4,  pp. 
343-353,  August  1978. 
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34.  J.  S.  Llm,  A.  V.  Oppenheim  and  L.  D.  Braida,  "Evaluation  of 
an  Adaptive  Comb  Filtering  Method  for  Enhancing  Speech  Degraded 
by  White  Noise  Addition,"  IEEE  Trans,  on  Acoustics,  Speech  and 
Signal  Processing,  vol.  ASSP-26,  no.  4,  pp.  354-358,  August  19 7 8". 

35.  J.  S.  Lim,  "Enhancement  and  Bandwidth  Compression  of  Noisy 
Speech  by  Estimation  of  Speech  and  its  Model  Parameters,"  Sc  ,D . 
Thesis ,  M.I.T.,  Dept,  of  Elec.  Eng.  and  Comp.  Science,  August 


36.  A.  V.  Oppenheim  and  J.  M.  Tribolet,  "Application  of 
Homomorphic  Filtering  to  Seismic  Data  Processing,"  chapter  in: 


Academic 


Time  Series  Analysis,  ed.  David  F.  Findley, 
Press,  pp.  26l-28(T,  September  1978. 


New 


York: 


37-  J.  S.  Lim,  "Evaluation  of  Correlation  Subtraction  Method  for 
Enhancing  Speech  Degraded  by  Additive  White  Noise,"  IEEE  Trans, 
on  Acoustics.  Speech  and  Signal  Processing  Letters,  vol.  ASSP-26, 
no.  5,  PP.  471-472,  "October  1978. 


38.  T.  F.  Quatleri,  Jr.,  "Short-Time  Spectral 
Conventional  and  Sliding  CZT,"  IEEE  Trans .  on 
and  Signal  Processing,  vol.  ASSP-26 , " no.  6 ,  pp . 
1978. 


Analysis  with  the 
Acoustics,  Speech 
561-566 ,  December 


1979 

39.  B.  R.  Musicus,  "An  Iterative  Technique  for  Maximum 
Likelihood  Estimation  with  Noisy  Data,"  S.M.  Thesis,  M.I.T., 
Dept,  of  Elec.  Eng.  and  Comp.  Science,  February  1979. 

40.  J.  S.  Lim,  "Spectral  Root  Homomorphic  Deconvolution  System," 
IEEE  Trans,  on  Acoustics,  Speech  and  Signal  Processing,  vol. 
ASSP-27,  no.  3,  pp.  223-233,  June  1979. 

41.  J.  S.  Lim  and  H.  Nawab ,  "Restoration  of  Speckle  Images," 
Lincoln  Lab.  Technical  Note  No.  1979-52,  July  3,  1979- 

42.  T.  F.  Quatleri,  "Minimum  and  Mixed  Phase  Speech 
Analysis-Synthesis  by  Adaptive  Homomorphic  Deconvolution,"  IEF.E 
Trans,  on  Acoustics,  Speech  and  Signal  Processing,  vol.  ASSP-27, 
no.  4,  pp. "328-335,  August  1979* 

43.  T.  F.  Quatleri,  Jr.,  "Phase  Estimation  with  Application  to 
Speech  Analysis-Synthesis,"  Technical  Report  No.  491,  Res.  Lab. 
of  Elec.,  November  1979. 

44.  J.  S.  Lim  and  A.  V.  Oppenheim,  "Enhancement  and  Bandwidth 
Compression  of  Noisy  Speech,"  IEEE  Proceedings,  (invited),  vol. 
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67,  no.  12,  pp.  1586-1604,  December  1979  (also  Lincoln  Lab. 
Technical  Note  1979-51.  June  20,  1979). 


1980 

45.  A.  V.  Oppenheim,  G.  V.  Frisk  and  D.  R.  Martinez,  "A 
Technique  for  the  Evaluation  of  Circularly  Symmetric 
Two-Dimensional  Fourier  Transforms  and  its  Application  to  the 
Measurement  of  Ocean  Bottom  Reflection  Coefficients,"  section  of 
a  book  entitled  DIGITAL  SIGNAL  PROCESSING,  eds.  V.  Cappellini  and 
A.G.  Constantinides ,  pp.  87-95,  London:  Academic  Press,  1980. 

46.  W.  P.  Dove,  "Event  Compression  Using  Recursive  Least  Squares 
Signal  Processing,"  Technical  Report  No.  492,  Res.  Lab.  of  Elec.. 
July  1980. 

47.  A.  V.  Oppenheim,  G.  V.  Frisk  and  D.  R.  Martinez, 
"Computation  of  the  Hankel  Transform  Using  Projections,"  Journal 
of  the  Acoustical  Society  of  America,  vol.  68.  no.  2.  pp. 
523-529,  August  1980. 

48.  G.  V.  Frisk,  A.  V.  Oppenheim  and  D.  R.  Martinez,  "A 
Technique  for  Measuring  the  Plane  Wave  Reflection  Coefficient  of 
the  Ocean  Bottom,"  Journal  of  the  Acoustical  Society  of  America, 
vol.  68,  no.  2,  pp.  602-612,  August  I98O. 

1981 

49.  A.  V.  Oppenheim  and  J.  S.  Lim,  "The  Importance  of  Phase  in 
Signals,"  IEEE  Proceedings  Special  Issue  on  Digital  Image 
Processing,  invited,  vol.  69,  no.  5,  pp.  529-541,  May  1981 . 

50.  J.  S.  Lim  and  H.  Nawab,  "Techniques  for  Speckle  Noise 
Removal,"  Journal  of  Optical  Engineering,  vol.  20,  no.  3.  PP. 
472-480,  May/June  19FT. 

51.  J.  S.  Lim  and  N.  A.  Malik,  "A  New  Algorithm  for 
Two-Dimensional  Maximum  Entropy  Power  Spectrum  Estimation,"  IEEE 
Trans .  on  Acoustics,  Speech  and  Signal  Processing,  vol.  ASSP-29, 
no.  3.  PP*  401-413,  June  1981,  (also  Lincoln  Lab.  Technical  Note 
1980-37,  Aug.  7,  1980). 
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PAPERS  SUBMITTED  FOR  PUBLICATION 

52.  M.  Hall,  A.  V.  Oppenheira  and  A.  S.  Willsky,  "Time-Varying 
Parametric  Modeling  of  Speech,"  IEEE  Trans,  on  Acoustics,  Speech 
and  Signal  Processing,  submitted  October  1979.  revised  September 
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